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We have calculated the fermionic contributions to the flavour non-singlet 
structure functions in deep-inelastic scattering (DIS) at third order of 
massless perturbative QCD. We discuss their implications for the thresh- 
old resummation at the next-to-next-to-leading logarithmic accuracy. 

1 Introduction 

The calculation of perturbative QCD corrections for deep-inelastic structure 
functions is an important task. The present and expected future experimen- 
tal precision, for instance at HERA, calls for complete next-to-next-to-leading 
order (NNLO) predictions. These offer the possibility to determine the strong 
coupling constant as and to analyze the proton structure and its parton content 
with unprecedented precision. Knowledge of the latter is of particular impor- 
tance for the analysis of hard scattering reactions at future LHC experiments. 

At present, this level of accuracy is not yet possible, because the necessary 
anomalous dimensions governing the parton evolution at NNLO are not fully 
known. The two- loop coefficient functions of F2 , F3 and F^ have been calcu- 
lated some time ago [1, 2, 3, 4, 5], but for the corresponding three-loop anoma- 
lous dimensions 7pp;, only a finite number of fixed Mellin moments [6, 7, 8] 
are available. As a first step towards the complete calculation, we have com- 
puted the fermionic three-loop contributions to the flavour non-singlet (NS) 
structure functions F2 and Fl in unpolarized electromagnetic deep-inelastic 
scattering [9]. Already these results have immediate consequences for thresh- 
old resummation of soft gluons which we will discuss in the following. 
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2 Threshold resummation 



It is well known that perturbative QCD corrections to structure functions re- 
ceive large logarithmic corrections, which originate from the emission of soft 
gluons. These corrections are relevant at large values of the scaling variable x 
(in Mellin space at large values of the Mellin moment N) and can be resummed 
to all orders in perturbation theory. It is interesting to investigate the implica- 
tions of our three-loop results [9] for the threshold exponentiation [10, 11, 12] 
at next-to-next-leading logarithmic (NNLL) accuracy [13]. 

Here the quark coefficient function for F2 can, up to terms which vanish for 
A'' — > 00, be written as 



C2{a,,N) = (l-hasffoi+«s5o2 + ---) exp[G^(Q2)], 
where the resummation exponent can be expanded as 

G^{Q'^) = L g^{as L) + g2{as L) + tts g^itts L) + . . . 
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with Us = as/(47r) and L = InA'^. The functions depend on universal coef- 
ficients Ai<i and Bi<i-i and process-dependent parameters Df^i_^ (see e.g. 
ref. [13] for the precise definitions of the functions g^ 2 3)- 

In a physical picture, the resummation of the perturbative expansion for 
C 2 rests upon the refactorization of C2 (valid in the threshold region of phase 
space) into separate functions of the jet-like, soft, and off-shell quanta that 
contribute to its quantum corrections. Each of the functions organizes large 
corrections corresponding to a particular region of phase space [14]. 

To NNLL accuracy, the function 93 involves the new coefficients A3, B2 

and D^^^. These coefficients can be fixed by expanding eq. (1) in powers of 

as and comparing to the result of the full fixed-order calculation [9]. In the 

MS scheme, the parameter A3 is simply the coefficient of In A'' in ^^ns{N) or, 

(2) 

equivalently, of 1/(1 — x)+ in Pna (x). It reads 
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where Ca = 3, Cp = 4/3 and is the number of hght (massless) flavours. 
The estimate of the non-fermionic part [13] is based on the approximations of 
Pns {x) constructed in ref. [15] using the first six even-integer moments [8] and 
its small-a; limit [16]. The exact fermionic part has been obtained independently 
in refs. [9, 18], while the nj^ contribution is ahcady known from ref. [17]. 

The complete results for B2 and D^^^ can be inferred from fermionic result 
of the three- loop coefficient function c 2"^^^ in ref. [9] , yielding 
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As a matter of fact, the contribution to c^^l^^ involves only a linear combination, 
/3o(^2 + 2 D^^^), with /?o being the one-loop coefficient of the QCD /?-function. 
However, the different combination B2 + D^^^ has been determined in ref. [13] 
by comparing the expansion of eq. (1) to the two- loop coefficient function c^'lys 
of ref. [19]. Thus, B2 and -Dj^'^ can be disentangled. It is interesting to observe 
the vanishing of D^^^ and -02^^^, for which an all-order generalization has been 
proposed in ref. [20, 21]. This is in contrast to the Drell-Yan process, where 
the functions are generally different from zero. For instance, has 

been derived in ref. [13]. 

Let us briefly illustrate numerically for large N the improvement due to the 
NNLL corrections for the soft gluon exponent of deep-inelastic scattering. 
In flg. 1 on the left, we show the resummation exponent {Q^) of eq. (2). 
Here, we choose fx^ = fj^f = Q"^, = A and as = 0.2, which corresponds 
to scales between about 25 and 50 GoV^, depending on the precise value of 
as(M|). In fig. 1 on the right, we display the convolution with a schematic, but 
typical input evaluated with the so-called 'minimal-prescription' contour [22]. 
It is obvious from both figures that knowledge of the leading logarithmic (LL) 
and next-to-leading logarithmic (NLL) terms [10, 11, 12] alone, i.e., those en- 
hanced by factors In A'^(q:s In A^)" and (aglnA^)", is not sufficient for reliably 
determining the function G^ and its impact after convolution even for rather 
large values of N and x. 

The NNLL corrections discussed here are rather small over a wide range. 
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This indicates that the soft-gluon exponent {Q"^) stabihzcs and that the soft- 
gluon effects on the MS quark coefficient function can be rehably estimated. 
Recall, however, that the NNLL corrections are large for the 'physical kernel' 
governing the scaling violations of the non-singlet structure function -F2,ns [23]. 
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Figure 1: Left: The LL, NLL and NNLL approximations for the resummation 
exponent {Q"^) in eq. (2) at /x^ = ^2 ^ q2 f^j. q;^(q2^ ^ o.2 and four 

flavours. Right: These results convoluted with a typical input shape xf = 
a;V2(i_^)3. 
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